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Abstract: Hereditary Tyrosinemia type 1 (HT1) is a metabolic liver disease caused by 
genetic defects of fumarylacetoacetate hydrolase (FAH), an enzyme necessary to complete 
the breakdown of tyrosine. The severe hepatic dysfunction caused by the lack of this 
enzyme is prevented by the therapeutic use of NTBC (2-[2-nitro-4-(trifluoromethyl)benzoyl] 
cyclohexane-l,3-dione). However despite the treatment, chronic hepatopathy and development 
of hepatocellular carcinoma (HCC) are still observed in some HT1 patients. Growing 
evidence show the important role of heat shock proteins (HSPs) in many cellular processes 
and their involvement in pathological diseases including cancer. Their survival-promoting 
effect by modulation of the apoptotic machinery is often correlated with poor prognosis 
and resistance to therapy in a number of cancers. Here, we sought to gain insight into the 
pathophysiological mechanisms associated with liver dysfunction and tumor development 
in a murine model of HT1. Differential gene expression patterns in livers of mice 
under HT1 stress, induced by drug retrieval, have shown deregulation of stress and 
cell death resistance genes. Among them, genes coding for HSPB and HSPA members, 
and for anti-apoptotic BCL-2 related mitochondrial proteins were associated with the 
hepatocarcinogenetic process. Our data highlight the variation of stress pathways related to 
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HT1 hepatocarcinogenesis suggesting the role of HSPs in rendering tyrosinemia-affected 
liver susceptible to the development of HCC. 

Keywords: Hereditary Tyrosinemia type 1 (HT1); fumarylacetoacetate hydrolase (FAH); 
NTBC (2-[2-nitro-4-(trifluoromethyl)benzoyl]cyclohexane-l,3-dione); hepatocellular carcinoma 
(HCC); heat shock proteins (HSPs); apoptosis 



1. Introduction 

Molecular chaperones play a key role as regulators of the apoptotic process [1]. Clinical data have 
demonstrated a correlation between increased HSPs expression and invasive potential of tumors [2,3]. 
Cancer cells are morphologically and functionally heterogeneous, their distinctive features make them 
able to survive despite an extreme environmental stress such as hypoxia, acidosis and elevated 
concentration of toxic metabolites. These stresses tend to generate free radicals that can damage 
cellular proteins. Since molecular chaperones have a protective role toward damaged proteins, it is not 
surprising to find them to be highly expressed in tumor cells. Most HSPs can behave as molecular 
chaperones for other cellular proteins participating in many cellular events with strong cytoprotective 
effects. HSPs have been classified into five families according to their molecular size: HSPH 
(HSP110), HSPC (HSP90), HSPA (HSP70), HSPD1 (HSP60), DNAJ (HSP40) and HSPB (small 
HSPs) [4]. Some HSPs are expressed under normal conditions (referred to as "constitutive" or 
"cognate") while others are induced during stress. Among molecular chaperones, HSPA and HSPB 
members are strongly induced by different stresses such as heat, irradiation, oxidative stress and 
anticancer therapy [5]. During carcinogenesis, the transformed cells begin to express an elevated level 
of HSPs and this induction continues during tumor progression [2,3]. At this time, it is still unclear if 
the over-expression of HSPs in cancer plays a role only in supporting malignancy or if it is essential in 
developing the transformed phenotype [1]. In the present study, we investigated the different 
expression patterns of HSPs and other anti-apoptotic proteins in liver cell transformation during 
hepatocarcinogenesis in a murine model of hereditary tyrosinemia type 1 (HT1). 

The liver plays a central role in the pathophysiology of many inborn errors of metabolism being the 
main site for catabolic, synthetic and detoxification reactions. Impaired degradation of the amino acid 
tyrosine is a feature of several genetic diseases mainly affecting the liver. Among them HT1 
(OMIM 276700) presents the most severe symptoms. HT1 is an autosomal recessive disorder 
characterized by severe liver damage, impaired coagulation, neurological crises, renal tubular 
dysfunction and a high risk of hepatocellular carcinoma [6-11]. The primary enzymatic defect in 
HT1 is a deficiency in fumarylacetoacetate hydrolase (FAH) activity, the last enzyme of tyrosine 
catabolism [12-17]. The lack of FAH produces an accumulation of the toxic upstream metabolites 
fumarylacetoacetate (FAA), maleylacetoacetate (MA A), and succinylacetone (SA) (Figure 1). 
This condition is responsible for the progressive injury to hepatocytes leading to chromosomal 
instability, and cell death in animals, and cultured cell models of the disease [18-23]. The only 
available treatment for this disease is the combination of a diet low in tyrosine and phenylalanine and a 
daily intake of NTBC (2-[2-nitro-4-(trifluoromethyl)benzoyl]cyclohexane-l,3-dione) also known as 
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nitisinone (Orfadin®, Swedish Orphan Biovitrum, Stockholm, Sweden). This drug inhibits 
4-hydroxyphenylpyruvate dioxygenase (HPPD), one of the enzymes involved in the catabolic pathway 
upstream of FAH, and therefore prevents the formation of toxic products responsible for liver damage 
(Figure 1). The efficacy of this therapy has improved the liver disease associated with FAH deficiency. 
However clinical data indicate that NTBC provides only partial protection against liver dysfunction. 
Indeed despite the improvement in survival and quality of life with NTBC treatment, HT1 remains a 
chronic disorder with several long-term complications, a persistent even though lower risk of HCC, 
and possible subnormal intellectual functions [24-27]. 

Figure I. The five enzymatic steps of tyrosine catabolism and NTBC site of action. 
NTBC acts as a pharmacological inhibitor of 4-hydroxyphenylpyruvate dioxygenase 
(HPPD) and is used in therapy to avoid production of homogentisate. TAT, tyrosine 
aminotransferase; HGO, homogentisate dioxygenase; MAAI, maleylacetoacetate isomerase; 
FAH, fumarylacetocetate hydrolase. 
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In mice HT1 is lethal but can be rescued by blocking upstream steps in the tyrosine catabolic 
pathway (i.e., NTBC). It has been reported that a chronic HT1 stress, induced by withdrawal of NTBC, 
results in resistance to cell death in vivo as shown by the survival of fah~ f ~ mice to an otherwise lethal 
dose of homogentisate [28]. This resistance led us to investigate how the balance between cell death 
and survival signals is maintained or disrupted in a murine FAH-deficient model of HT1. Therefore, 
we analyzed the relationship between the evolution of major survival and apoptotic pathways and the 
progression of clinicopathologic features of the disease. Our previous studies showed that the HT1 
phenotype induced by NTBC-withdrawal causes a cellular insult eliciting the ER stress response and 
increasing the level of stress-related proteins [29]. This hepatic stress causes the activation of many 
survival pathways and inhibits the intrinsic apoptotic cascade promoting hepatocarcinogenesis [30]. 
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Thus, to clarify the biologic pathways that degenerate the HT1 stress response to cancer, we examined 
the gene expression patterns in tumoral liver fractions compared to healthy livers. Our studies show 
how the survival state in HT1 livers under progression of cancer promoted by long-term NTBC 
withdrawal, can influence the expression of HSPs and other pro-survival proteins that act at the 
mitochondrion level. Indeed, molecular chaperones have multiple roles in cell survival, that depend on 
their distinctive features, one of them being their interactions with anti-apoptotic proteins of the BAG 
and BCL-2 families [31]. 

2. Results and Discussion 

2.7. Gene Expression Analysis infaK f ~ Mice after NTBC Withdrawal 

The murine FAH-deficient model of HT1 presents a neonatal lethal phenotype that can be rescued 
by treating the pregnant and nursing females and the pups with NTBC [32]. FAH deficient mice have a 
severe phenotype, displaying rapidly a remarkable impairment of vital functions upon withdrawal of 
the drug. Tht fah~ y ~ mice receive NTBC in drinking water at 7.5 mg/L from birth. Assuming that they 
drink 3-5 mL of water per day, the mice are kept on a high-dose treatment. Our previous studies 
demonstrate the changes in the activation state of cellular signaling pathways in response to HT1 stress 
that have been observed after one to five weeks of NTBC withdrawal [30]. With the aim to assess the 
different patterns of genes expression during a long-term HT1 stress characterized by tumor development, 
we performed differential gene expression profile in livers of tyrosinemic mice having developed HCC 
following cumulative HT1 stress after a period of NTBC interruption from 5 to 7 weeks relative to 
livers of wild type and/a/T /_ control mice (Figure 2). Selective dissection of tumor and non-tumor 
tissues and extraction of total RNA allowed the analysis of gene expression with DNA microarrays of 
Whole Mouse Genome (Agilent Technologies, Santa Clara, CA, USA). We investigated the 
modulation of genes and their change in tumoral fractions compared with fah +/+ and NTBC-treated 
fah~ f ~ mice. Differently expressed genes were selected by filtering with confidence at p < 0.05, 
considering those with a difference in level of expression of at least 2-fold. Biological classification of 
differently regulated genes was obtained by Gene Ontology (GO) analysis (PANTHER classification 
system). Tumoral samples displayed significant modulation in the expression of genes leading to cell 
survival. The GeneSpring algorithm (Agilent) and the IPA software (Ingenuity, QIAGEN, Redwood 
City, CA, USA) were further used to create functional classes of genes implicated in regulation of the 
survival pathway in the context of HT1 (Table SI, Supplementary Data). The wide range of genes 
deregulated in the tumoral fraction included a substantial number of genes implicated in cellular 
growth and proliferation, organismal survival and apoptosis, differentiation, inflammation and cell 
migration (Figure 3 and Table SI, Supplementary Data). These results show that in tumor cells a 
preferential expression of specific constellation of genes related to hepatocarcinogenesis occurs. 

2.2. HSPs Expression Is Associated to HT1 Progression in Mice 

Notably, among genes deregulated in tumoral fractions we found those belonging to HSP families. 
Liver samples under HT1 stress showed over-expression of HSPs transcripts; among them, Hspbl and 
the HSP70 member Hspala and Hspa2 were the most up-regulated compared to wild type (Figure 4A) 
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and control mice (Figure 4B). Following NTBC withdrawal, Hspala showed a 4.61 fold increase in 
fah~'~ NTBC-w/o mice compared to fah +/+ mice and 3.28 fold compared to fah / NTBC-treated mice. 
Hspbl was also higher in mice after NTBC withdrawal presenting fold changes of 2.36 and 2.77 
compared to fah +/+ and NTBC-treated fah~ f ~ respectively (Figure 4). Other transcripts of the HSP 
family were also induced such as Hspa2, Hspb3 and DnajclO. Hspa2 encodes heat shock-related 70 kDa 
protein 2 involved in normal development of germ cell and is over expressed in cancer 
cells [33,34], bladder urothelial carcinoma [35] and peptic ulcer [36]. 

Figure 2. (A) Biological processes regulated by all significant differentially expressed 
genes assessed by Gene Ontology (GO) search and summarized according to their 
functions (PANTHER classification system); (B) The GO level 1 of the response to 
stimulus was further mined down showing as subset the percentage of molecules involved 
in stress response. 
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Figure 3. Disease and function annotation according to IPA software. As it can be noticed 
the largest amount of deregulated genes is involved in the mechanism of cancer. Others 
relevant processes are those implicated in gastrointestinal disease, cellular growth and 
proliferation, migration and differentiation of cells. These evidence suggest that under HT1 
stress there is a wide gene regulation involving many biological processes that all together 
might foster the tumoral progression. 
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Figure 4. HSPs genetic expression. Influence of HT1 stress on expression levels of 
transcripts codifying for HSP70 family (Hspala, Hspa2), HSPB family (Hspbl, Hspb3) 
and DNAJ co-chaperone family. Graphics showing the relative expression level in livers 
after withdrawal of NTBC in/a/T^mice compared to livers of wild type (A), and of fah~ f ~ 
NTBC-treated mice (B). Data are given as mean of two experiments (n = 2) in a single 
scan. Each sample is representative of a pool of four different mice for each condition. For 
all up-regulated genes p < 0.05 was obtained by variance analysis (ANOVA). The tumoral 
fractions showed an up regulation in the HSPs families up to 5-fold higher compared to the 
fah +/+ (A) mdfah'^ control (B) mice. 
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The Hspb3 gene encodes for HSPB3 (HSPL27), a small heat shock protein involved in distal motor 
neuropathy [37]. Unlike the other HSPs members, Hspb3 was significantly deregulated in fah~ f ~ 
tyrosinemic mice compared to fah +/+ wild type mice, but no difference was observed when compared 
to fah~ f ~ NTBC-treated mice (Figure 4). DnajclO is a member of the HSP40 family and encodes for 
the endoplasmic reticulum resident protein DNACJ10 shown to be up-regulated after ER-stress [38,39]. 
Since the functions of HSPs in tumorigenesis have been related to their ability to interfere with the 
apoptotic machinery, we also examined the level of mRNAs encoding anti-apoptotic proteins that act 
at the mitochondrial level. A subset of transcripts of the BCL-2 family was found to be over-expressed in 
tumors of HT1 mice (Figure 5). In particular the anti-apoptotic members of the BCL-2 like family, 
Bcl2ala, Bcl2alb, Bcl2lll and Bcl2, showed a remarkable fold change compared to fah +/+ and /a/T /_ 
NTBC-treated (Figure 5). Moreover, Bag2 gene, member of the BCL-2 athanogene related proteins, 
showed an expression value respectively of 5.49 fold and 4.73 fold higher in wild type (A) and 
NTBC-treated mice (B). 

Figure 5. BCL-2 related protein mRNAs expression in tumoral liver tissue. Comparison of 
gene expression values of the BCL-2 family members between healthy and tumoral livers. 
Graphics showing the relative expression level in livers after withdrawal of NTBC in 
fah~ f ~ mice compared to liver of wild type fah +/+ mice (A), and/a/T /_ NTBC-treated mice 
(B). Data are given as mean of two experiments (n = 2) in a single scan. Each sample is 
representative of a pool of four different mice for each condition. For all deregulated genes 
p < 0.05 was obtained by variance analysis (ANOVA). A ratio up to 9-fold higher of the 
BCL-2 related protein mRNAs is observed. In particular the anti-apoptotic Bcl2ala, Bcl2alb, 
Bcl2, Bcl2lll, were over-expressed, while the pro-apoptotic Bmf was down-regulated. 
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2.3. Increase of Pathological Feature during Long-Term HT1 Stress Promotes HCC Development 

In the case of the NTBC withdrawal induced HT1 stress, the physiological state of the/a/T 7- mouse 
reflects a constant disease progression, proportional to the damage accumulation. Several studies on 
the role of HSPs in tumorigenesis have established that inducible HSPs are up-regulated in cancer 
cells. HSPs have also been involved in facilitating cell survival through the inhibition of 
apoptosis [40], and cell senescence [41]. Hence the result of gene expression analysis prompted us to 
investigate the participation of molecular chaperones in the biological processes that degenerate the 
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HT1 stress into cancer. Thus, we built a new protocol with the purpose of inducing a progressive 
severe phenotype of the disease in which 4-months old fah~ ; ~ mice were withdrawn from NTBC 
treatment for periods up to 15 weeks. Figure 6 shows the time points of the protocol chosen to evaluate 
the different grades of injury post-withdrawal. 

Figure 6. Protocol of study. 
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Four time points were set to look at the disease in different stages until formation of large tumors in 



100% of the mice at 15 weeks. Young fa\f f ~ mice were treated from birth until the age of four months 
with NTBC. As control, fah +/+ healthy mice were used. After four months of NTBC treatment the drug 
was removed to allow HT1 development and we started the harvest. The mice were examined daily 
and any variation of physiological conditions was recorded. The animals were weighed three times a 
week until the end of the protocol (15 weeks). As previously reported, the fah knockout mice 
displayed a severe weight loss and a rapid hepatic dysfunction following the drug withdrawal. After 
five weeks without NTBC, fah~ f ~ mice display a loss of total body mass reaching -38% from initial 
mass (Figure 7A). During the next 10 weeks they regain weigh up to -80% of the original body mass 
until the time of sacrifice. This gain of weight was mainly due to the liver that doubled its mass and 
presented many macronodules (Figure 7B,C). The resulting slope of total weigh was a sign of clinical 
illness leading to poor prognosis. Indeed, all mice harvested at the end of the protocol had the tendency 
to develop a severe tumor phenotype (Figure 7C). Histological analysis of the livers at the different 
time points showed a progressive hepatic pathology with inflammation, severe histological lesions, 
architectural rearrangements, severe hepatocellular changes, apoptosis, and many erythropoietic foci 
(data not shown). 

2.4. Variation in the Expression ofHSPAlA and HSPB1 Is Correlated to HT1 Progression in 
NTBC-withdrawn Mice 

Liver homogenates obtained from mice taken off NTBC from one to 15 weeks, were next used to 
determine the level of stress induction and resistance to apoptosis. HSPC2/3 (HSP90 a/(3), HSPD1 
(HSP60) and HSPA5 (BiP) showed a constant level of expression during disease progression 
(Figure 8A). The cognate form of the 70kD heat shock protein, HSPA8, remains almost unvaried in 
fah~ f ~ mice, independently of the NTBC treatment, with a faint reinforcement in the advanced 
hepatocarcinogenetic stage. In contrast, the inducible HSP70 member HSPA1A and the small HSPB1 
were significantly deregulated along the entire experience (Figure 8A,B). 
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Figure 7. Changes in body mass values (%) of wt (A) and fah ' (■) mice following 
NTBC withdrawal. (A) Body weight curves of fah +/+ (A), and /a/T /_ NTBC w/o (■) mice. 
Each time point indicates the average of four mice, with error bars indicating standard 
deviation. Unpaired t-test (Welch's correction) was applied for statistical validation of the 
two groups (p value < 0.0001 for all the time point). (B) Liver weight (a), and liver /body 
weight ratios (b), of fah +/+ , fah~ f ~ NTBC-treated mice andyh/T 7- mice following 15 weeks 
of withdrawal (n = 4, ****^< 0.0001, ANOVA variance analysis). (C) Representative 
pictures of whole healthy livers and tumoral livers. The images were taken at the beginning 
of the procedure (control mice) and after 15 weeks of withdrawal from treatment. 
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Figure 8. Western blot analysis of heat shock proteins. (A) Variation in the levels of 
anti-apoptotic HSPB1 and HSPA1A, in liver of fah +/+ , fah'^ NTBC-treated and fah'^ 
NTBC-off mice. HSPC2/3, HSPA5, HSPA8 and HSPD1 were also tested as HSP family 
members. Data are representative of at least two different experiments with four biological 
replicates; (B) Histograms show quantitative representation of the level of HSPs proteins 
after normalization to levels of GAPDH and the fold-change protein level expression is 
reported in comparison to fah +/+ mice. Data are mean + standard deviation for four 
samples. Significance of group differences was evaluated using ANOVA for multiple 
comparisons (*** p < 0.001 and **** p < 0.0001). Where statistical significance was not 
observed after densitometric analysis, graphics were not produced (i.e., HSPC2/3, HSPA5, 
HSPA8, HSPD1). 
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HSPA1A shows a remarkable variation in its expression during the starvation from the drug. It is 
already expressed in fah +/+ and fah~ f ~ NTBC-treated mice. Soon after the withdrawal it goes up for the 
first week, then it starts to come back to the original level (Figure 8A,B(a)). HSPB1 and HSPA1A 
proteins have an anti-apoptotic role in cells by inhibiting the formation of the apoptosome, 
thus eliciting the intrinsic apoptotic pathway. Induction of HSPB1 appears already three days after 
NTBC-treatment interruption and increases notably during HT1 stress progression with a 
reinforcement at four weeks of withdrawal. Its phosphorylated form on serine 15 appears also 
increased in the same proportion (Figure 8A,B(b,c)) which is different from what has been reported in 
some experimental studies where an attenuated phosphorylation of HSPB1 in advanced HCC was 
observed [42,43]. The phosphorylated form of HSPB1 at serine 15 is involved in the ability of this 
chaperone to switch between small and large oligomers in order to accomplish its anti-apoptotic 
functions [40,44-47]. As reported by other authors the equilibrium between large and small oligomers 
might be shifted by phosphorylation, depending on the physiological requirements of the cell. 
Interestingly, the ability to form small oligomers upon phosphorylation was shown to mediate the 
entry of HspBl in the nucleus [48-51]. Based on these findings HSPB1 phospho-Ser-15 might 
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perform other important functions in HT1 -dependent tumorigenesis by entering into the nucleus as 
small oligomer. Further studies will be needed to investigate the crosstalk between aggregated and 
small-dissociated forms that could regulate HSPB 1 sub-cellular localization and biological functions in 
liver cancer during HT1 stress. In this murine model the chaperone-induced cytoprotection could 
rescue cells from apoptosis, easing the deleterious consequences of HT1 chronic stress. It is also 
important to notice the variability of expression in the HSPs family in the different stages of the 
disease. In fact, as previously reported, HSP70 and HSPB 1 could be involved in different phases of 
carcinogenesis implicating a complex network of apoptotic pathways, supporting the survival state and 
favoring the progression of the malignancy. Tumoral cells could use the cytoprotective effect of 
HSP70 and HSPB1 and their ability to interact with many targets, to survive to the HT1 -driven 
environmental stress promoting the tumoral invasion. 

2.5. Survival State in the Liver Is Reinforced by Increased Expression of the Anti-Apoptotic 
BCL-2 Protein 

In agreement with the microarray result and with our previous studies [30], the levels of the 
anti-apoptotic BCL-2 protein were also affected (Figure 9). BCL-2 is the central coordinator of 
apoptotic events into the mitochondria. The rupture of the outer mitochondrial membrane induces the 
mitochondrial translocations and multimerization of the pro-apoptotic BCL-2-associated X protein 
(BAX) resulting in the leakage of cytochrome c and other mitochondrial proteins for the execution 
phase of apoptosis. In this context BCL-2 plays a fundamental role by preventing the mitochondrial 
generation of ROS and dimerizing with BAX, thus blocking the BAX interaction with the 
voltage-dependent anion channel (VDAC) to form pores that are large enough to allow the passage of 
cytochrome c [52,53]. In the HT1 murine model the protein levels of the anti-apoptotic BCL-2 show a 
remarkable variation with the progression of the disease (Figure 9). Indeed, BCL-2 was undetectable 
in fah +/+ mice, but it is already expressed in the fah~ f ~ control, increasing its expression in the first 
week up to four weeks without treatment, and then starting to decrease its level again (Figure 9A,B). 
The trend of BCL-2 seems to be comparable to HSPA1 A as it tends to increase with the progression of 
the disease but, in contrast to HSPA1A, it seems to continue until the fourth week of withdrawal and 
then decrease to the beginning level. Over expression of BCL-2 has already been found in various 
types of cancer [54] and it is considered as a biomarker of resistance to chemotherapy and 
radiotherapy [55-57] although its universal role in cancer is still unclear [58]. These results suggest the 
involvement of BCL-2 together with the HSPs, in the biological series of events that exacerbate the 
HT1 stress leading to hepatocarcinogenesis. Herein we could conclude that in the first phase of HT1 
stress BCL-2 contributes, cooperating with HSPs, in the apoptosis resistance that facilitates tumor onset. 

2.6. HSPs and Tumor: An Intricate Co-Operation 

From these results we suggest that HT1 related liver cancer could arise from a complex network of 
gene deregulation due to the metabolic defect. Activation of many survival pathways has already been 
reported in HT1 cells models [20], and resistance to cell death was also documented in HT1 
stress [28,29]. In these studies an induction of HSPB1 and HSPA1A protein levels has been observed 
mfah~ H mice at different periods after NTBC withdrawal. It has been shown that HSPs can participate 
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in multiple cellular processes pertinent to tumor aggressiveness promoting cancer development by 
increasing cellular migration [59] and differentiation [60]. Although these HSPs are not mediators of 
proliferation, their role in tumorigenesis involves their ability to block programmed cell death (PCD). 
They have been shown to play an essential role in the survival of a wide number of human cancer cells 
by interfering with multiple apoptotic pathways. HSPA1A and HSPB1 have a strong cytoprotective 
effect that increases metastatic potential and growth ability in tumor cells, as a consequence of the 
inhibition of the apoptotic machinery [61,62]. These HSPs can regulate in an independent manner the 
apoptotic cascade upstream of the mitochondria by inhibiting stress inducible signals, at mitochondrial 
level by preventing membrane permeabilization and release of cytochrome c, and by suppressing 
caspase activation downstream the mitochondria [63]. HSP70 has also been implicated in other 
mechanisms of cell death in addition to PCD. This HSP has been shown to increase the lysosomal 
resistance against chemical and physical membrane destabilization [64]. Nevertheless it has been 
suggested that HSP70 could modulate JNK and ERK phosphorylation in apoptosis induced by 
hyperosmolarity [65] and promote the stabilization of phosphorylated form of PKC [66]. 

Figure 9. BCL-2 protein expression in HT1. (A) Levels of BCL-2 in livers of fah +/+ ,fah _/ ~ 
NTBC-treated mdfah'^ NTBC-off mice; (B) Histogram shows quantitative representation 
of the BCL-2 protein level after normalization to levels of GAPDH. Fold-change protein 
level expression is reported in comparison to fah +/+ mice. Data are mean + standard 
deviation for four samples. Significance of group differences was evaluated using ANOVA 
for multiple comparisons (** p < 0.01, *** p < 0.001, **** p < 0.0001). 
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HSPB1 is an ATP-independent chaperone and interacts with its substrates by forming multimeric 
complexes [67]. The oligomerization of HSPB1 is a dynamic process highly regulated by 
phosphorylation at different serine residues (Serl5, Ser78, Ser82 in humans; Serl5 and Ser86 in mice) 
and involved in the determination of substrate affinity. Experimental evidence suggest that HSPB1 
mainly acts as an inhibitor of caspase-dependent apoptosis and this action seems to be related with 
large unphosphorylated oligomers [45,47]. However the direct interaction with cytochrome c does not 
require the phosphorylation of HSPB1 [44]. Still, small-phosphorylated oligomers of HSPB1 
efficiently bind Daxx and F-actin enhancing its cytoprotective functions [40,46]. 
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Consistent with all these observations, these inducible chaperones have been implicated in 
cancer progression and in chemotherapy resistance and their overexpression is considered 
as negative prognostic marker for different cancers, as colorectal cancer [68-70], intrahepatic 
cholangiocarcinoma [71], prostate cancer [72] and urothelial carcinoma [73]. Furthermore, with the 
aim to find new strategies to block tumoral invasion, HSPB1 and HSP70 have been proposed as target 
for anticancer therapy [74,75]. Our results show that transcripts encoding specific HSPs and 
the corresponding proteins are induced in HT1 livers during carcinogenesis. This might lead to 
further investigations to understand the role that these HSPs may have as prognostic biomarkers 
for HCC development in HT1 patients and in developing new therapies in addition to NTBC treatment 
in these patients. 

HSPs and BCL-2 families regulate cell survival through many interactions that depends on their unique 
features. Tumor cells may escape from caspase-mediated apoptosis by overexpressing anti-apoptotic 
proteins. Thus our investigations were also extended to the level of proteins that control the cell's 
response to apoptosis. Regulation of the mitochondrial pathway of apoptosis via BCL-2 family 
members is supported by gene and protein expression. BCL-2 family consists of about 25 members 
and includes pro-apoptotic factors such as BAD, BAK, BAX, BID, BIM and anti-apoptotic proteins 
including BCL-2 and BCL-X L . The balance between members of BCL-2 family governs cell's response 
to apoptosis signaling. Altered expression of the survival genes such as those codifying for BCL-2 
anti-apoptotic proteins induces dysfunction in the apoptotic machinery and has therefore been associated 
with carcinogenesis [76-82]. Among genes up-regulated in HT1 tumoral fractions we noticed also one 
transcript of the BAG family member, BAG-2 (Figure 5). The BAG family (BCL-2 associated athanogene) 
represents nucleotide exchange factors in the ATP-dependent chaperone cycle of HSP70 [83-85]. They 
make a direct link between HSP and BCL-2 helping BCL-2 activation [86]. BAG family members 
possess an ubiquitin-like domain that allows their association with the proteasome [87,88]; this cooperation 
not only stimulates chaperone-assisted degradation but also interferes with it. In particular, BAG2 was 
reported to inhibit the CHIP-mediated ubiquitylation of HSP70, despite other BAG co-chaperone 
members, such as BAG-1, that promote the release of HSP70 in the vicinity of the proteasome 
stimulating the HSP70 degradation [89]. We need to investigate more about this involvement in HT1 
stress but it is nevertheless tempting to speculate that the cooperation between HSPs, BCL-2 and BAG 
family members exerts regulatory function contributing to the cell death resistance that promote 
carcinogenesis in HT1 stressed cells. 

3. Experimental 

3.1. Animal Maintenance and Treatments 

HT1 phenotype was reproduced using homozygous FAH deficient mice (fah AQXOn5 , referred as 
fah~ f ~) [27]. All animals were originally inbred of 50% 129SvJ and 50% C57BL/6 strains. The j ah 
knockout mice and normal mice fah +/+ were maintained in a temperature and humidity-controlled 
environment with a 12 h light-dark cycle. Mice were allowed access to rodent food (0.82% tyrosine, 
0.51% phenylalanine, Purina 5075-U.S., Charles River, St-Constant, QC, Canada) and drinking water, 
which was supplemented with 2-(2-nitro-4-trifluoromethylbenzoyl)-l,3-cyclohexanedione, NTBC 
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(kindly provided by S. Lindstedt, Gottenburg University, Sweden) for the fah~ ~ mice. The NTBC was 
used at a concentration of 7.5 mg/L and provided from fetal stage by treatment of the pregnant and 
nursing females. To induce the HT1 phenotype mice were withdrawn from NTBC treatment for 
different periods of time from 1 to 15 weeks. Mice control fah +/+ and mice fah~'~ on NTBC were used 
to underline the differences dues to genotypes. Mice were treated according to the guidelines of the 
Canadian Council of Animal Care (CCAC) and euthanized at the time of sacrifice by cardiac puncture 
with an overdose of ketamine/xylazine (0.2 mL/lOg). 

3.2. Gene Expression Analysis 

Total RNA was isolated from snap-frozen liver tissues using TRIzol Plus Reagent (Invitrogen 
Corporation, Burlington, ON, Canada) according to the manufacturer's instructions. Total RNA was 
treated with DNase using the RNase-free DNase kit and RNeasy spin column (QIAGEN, Toronto, ON, 
Canada) and dissolved in RNase-free water to a final concentration 0.2-0.5 |ig/|iL. The quality and 
purity of RNA samples was checked by capillary electrophoresis using the Agilent 2100 Bioanalyzer 
and Agilent RNA Nano 6000 LabChip kits (Agilent Technologies). Three different RNA pools 
resulting from four different mice livers were used to perform the hybridization on the four arrays. 
Five hundred nanograms of purified RNA were used to prepare cDNA labeled with Cy3 and Cy5, 
using the two colors Agilent Low Input Quick Amp Labeling kit, and hybridized using Agilent 
Two-Color Microarray-based gene expression platform using the 4x44K Two Colors Mouse Whole 
Genome Microarray, 41,534 genes (Agilent Technologies). Two hybridizations were carried out for 
tumor against reference samples using a fluorescent dye reversal technique (dye-swap). Slides were 
scanned using an Agilent DNA Microarray Scanner and Cy3/Cy5 intensity data were extracted using 
Agilent's Feature Extraction software (Agilent Technologies). Transcriptome analysis was performed 
using Agilent Genespring GX11 software (FRSQ-Reseau Cancer Core Genome Facility, Hotel-Dieu, 
Quebec, QC, Canada). A list of differently expressed genes was first created, by filtering on 
confidence at p < 0.05, consecutive lists were generated considering a fold change value more then 
2-fold. IPA Ingenuity software (QIAGEN, Redwood City, CA, USA), was further used to create 
functional classes of genes implicated in regulation of the survival pathway in the context of HT1. 
Biological processes regulated by all significant differentially expressed genes were classified 

TM 

according to their function using the Gene Ontology (PANTHER GO slim). 

3.3. Protein Electrophoresis and Western Blot Analysis 

Liver tissues were immediately snap frozen in liquid nitrogen and stored at -80 °C. For western 
blot, frozen tissues were homogenized in lx RIPA Buffer (20 mM Tris-HCl, pH 7.5; 150 mM NaCl; 
1 mM Na 2 EDTA; 1 mM EGTA; 1% NP-40) with IX protease inhibitor cocktail tablets (cOmplete, 
Mini, EDTA-free; Roche Diagnostics, Indianapolis, USA) and lx phosphatase inhibitor cocktail 
tablets (PhosSTOP; Roche Diagnostics, Indianapolis, IN, USA). The homogenized extracts were 
centrifuged at 10,000 g for 15 min at 4 °C. Protein concentration was measured with the Bio-Rad 
Protein Assay (Bio-Rad Laboratories Inc., Hercules, CA, USA). Proteins were resolved by 
electrophoresis on a 12% SDS-PAGE gel, transferred onto a nitrocellulose blotting membrane 
(BioTrace NT, Pall Life Ltd., St-Laurent, QC, Canada), blocked with 5% non-fat dried milk in TBS 
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(Tris-buffered saline; 50 mM Tris, 150 mM NaCl, pH 7.5) and probe overnight at 4 °C with primary 
antibodies following recommendations by the manufacturer. Home-made rabbit polyclonal antibodies 
against HSPA1A (#799 1:10,000), HSPA8 (#1477 1:10,000), HSPC2/3 (#214 1:10,000) and HSPD1 
(#37 1:10,000), were used and are described in Tanguay et al. [90]. Rabbit polyclonal antibody against 
HSPB1 was kindly provided by J. Landry (#L2R3 1:1000, Hotel-Dieu Research Center, CHU-Q, 
Quebec, QC, Canada). Rabbit polyclonal antibody against HSPB1 phospho-Serl5 and HSPA5/BiP 
were purchased from Stressgen Biotechnologies Corp. (Victoria, BC, Canada) (cat. #SPA-525 1:1000, 
and #SPA-826 1:1000). Monoclonal rabbit antibodies against BCL-2 and GAPDH (cat. #2870 1:1000 
and #5174 1:2000) and secondary anti-rabbit IgG-HRP (cat. #7074 1:2000) antibody were from Cell 
Signaling Inc. (Beverly, MA, USA). Odyssey® Infrared Imaging System (Li-COR, Biosciences, 
Lincoln, NE, USA) was used for proteins detection. Densitometric analysis were performed on 
western blot images using ImageJ 1.47v. The obtained values were normalized to GAPDH and finally 
were compared to the control fah +/+ mice. Data in histograms indicate the average of four mice for 
each time point, with error bars indicating standard deviation. 

3.4. Statistical Analysis 

GeneSpring GX11 software (Agilent Technologies) was used to generate a list of selected genes 
and for different statistical methods on gene expression data. The significance of group differences 
was evaluated using ANOVA multivariance analysis, followed by post hoc Tukey test for selection of 
genes in pairwise comparison. Unpaired Student's t-test was used in comparing weight variation of 
fah +/+ and fah~ f ~ NTBC- withdrawn mice during the protocol of study. Ordinary one-way ANOVA 
variance analysis (Brown-Forsythe test) was applied to measure the statistical significance of the trend 
of liver weight and liver/body weight ratio in the different groups. Statistical evaluation on western 
blot images was performed by two-way analysis of variance (ANOVA) followed by Dunnett's 
multiple comparison test (GraphPad Prism 6.0d). A statistically significant difference was defined as 
* p < 0.05, ** p < 0.01, *** p < 0.001, **** p < 0.0001. All p values were derived from two-tailed tests. 

4. Conclusions 

In summary, our results demonstrate that HT1 stress induces progressive damage on hepatic system 
as soon as the HT1 stress is initiated by NTBC withdrawal, concomitant to an early, constant and 
specific activation of the survival pathway in livers of fah~ f ~ mice. In early and in advanced 
carcinogenesis, these pathways seem to be differently modulated suggesting that more detailed 
investigations are needed to elucidate whether they are simply fast-acting proteins in response to tumor 
stress or indeed, specific marker for different phases of cancer progression and development. 
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